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Abstract. The present paper is devoted to the detailed study 
of quantization and evolution of the point limit of the Pauli-Fierz 
model for a charged oscillator interacting with the electromagnetic 
held in dipole approximation. In particular, a well defined dynam- 
ics is constructed for the classical model, which is subsequently 
quantized according to the Segal scheme. To this end, the classical 
model in the point limit, already obtained in |20|. is reformulated 
as a second order abstract wave equation, and a consistent quan- 
tum evolution is given. This allows a study of the behaviour of the 
survival and transition amplitudes for the process of decay of the 
excited states of the charged particle, and the emission of photons 
in the decay process. In particular, for the survival amplitude the 
exact time behaviour is found. This is completely determined by 
the resonances of the systems plus a tail term prevailing in the 
asymptotic, long time regime. Moreover, the survival amplitude 
cxhibites in a fairly clear way the Lamb shift correction to the 
unperturbed frequencies of the oscillator. 



1. Introduction 

In recent years a considerable effort was tributed by the mathemat- 
ical physics community to the problem of a rigorous formulation of 

l 
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the dynamics of the main models in nonrelativistic quantum field the- 
ory. In particular, a comprehensive study of the Pauli-Fierz model, 
the model which describes the low energy interaction of nonrelativis- 
tic matter and electromagnetic radiation, was undertaken by various 
authors, both in its full form, or making resort to different approxi- 
mations, such as rotating wave approximation, dipole approximation, 
or others. Correspondingly, a wealth of results concerning various as- 
pects of the model have been obtained, concerning self-adjointness of 
the hamiltonian, existence, multiplicity or also non existence of the 
ground state and related infrared behaviour, and detailed study of the 
spectral properties of the model and of its resonances. In the present 
paper we give a comprehensive analysis of some of the previous prob- 
lems in the special case of a point charged oscillator interacting with 
the electromagnetic field in dipole approximation. While this model 
is unrealistically simple compared to the case, to give an example, of 
the hydrogen atom interacting with the full (not dipole) radiation field, 
(about which a lot is known thanks to the work of Lieb et al. [T2].[T6] 
and Frolich et al. [H],[I], the novelty of the present work resides in the 
fact that we are able to cope with the point limit of the model. The 
removal of the ultraviolet cutoff in the interaction between matter and 
radiation is in its generality, a difficult and unsolved problem, and in 
particular one not faced off in the quoted rigorous literature. In some 
previous papers of the last two authors ([131201 and references therein) 
the renormalized dynamics of the classical Pauli-Fierz model in dipole 
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approximation ([211) an d f° r fairly general external potentials, was rig- 
orously constructed and analyzed. In particular, it was shown that the 
evolution of the Pauli-Fierz model in the point limit is given by an 
abstract wave equation generated by a family of operators related to 
the so called point interactions (see pQ and references therein). Given 
the classical model in the ultraviolet limit, a second step would be to 
construct the quantized model. The harmonic potential has the unique 
feature of giving rise to linear equations of motion, so that the classical 
equations for the system take the form of an abstract linear wave equa- 
tion. This allows a plain quantization a la Segal of the model, which 
seems otherwise quite problematic (in contrast with the case of the 
regularized Pauli-Fierz model, where canonical quantization works: a 
detailed study of the regularized dipole Pauli-Fierz model with an ex- 
ternal harmonic potential was given by Arai in [2 ). So we confine 
ourself to this case, about which we give in the next section a self- 
contained treatment independent and different in spirit from [20), to 
clarify some of the themes discussed above. We outline briefly our main 
results. The classical evolution of the system is given by the abstract 
wave equation 

(dl + L e )£ = 

for the couple £ = (A,p) where A is the vector potential of the elec- 
tromagnetic field in Coulomb gauge and p is the particle momentum 
variable (see section 2 for the explanation of this choice). Here L e is 
a self-adjoint operator in L^(IR 3 ) © R 3 (* stands for "divergenceless" ) 
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such that its resolvent can be explicitly calculated (see Lemma 2.2 and 
Theorem 2.3). 

The operator L e has an single negative eigenvalue and the rest of its 
spectrum is purely absolutely continuous and coincides with [0, +oo). 
Thus, to quantize the abstract wave equation above according to the 
Segal method, one has to take L+, i.e. L e projected onto the spectral 
subspace corresponding to cr ac (L e ). The (first) quantized dynamics 
of the system is defined through the Schrodinger-like equation iip = 
(L+y/ 2 ^, defined on the complex Hilbert space L^(R 3 ; C 3 )©C 3 and the 
hamiltonian of the quantized system of particle and field is given by the 
second quantization c£T((L+)5), on the Fock space over L*(R 3 ; C 3 )©C 3 . 

A preliminary but essential step in the description of the proper- 
ties of the system is to write the first quantized evolution in terms of 
the resolvent of the original classical operator L+. By Stone formula, 
spectral theorem and after some work, one gets (see Lemma 2.5) a 
representation (perhaps new or at least not known to us) for the tran- 
sition amplitudes (ipi, e~**( L ^) 1/s ''ip 2 ) between one particle states in terms 
of the resolvent of L e . Correspondingly, one has an expression for the 
survival and transition amplitudes for the second quantized model on 
the Fock space, just using functoriality of T. Our main concern here 
is in the calculation of two relevant characteristics of the evolution. 
The first is the survival amplitude S (t) of the unperturbed first excited 
bound state of the oscillator. The second is the amplitude transition 
A(t) between states with the photon field in the vacuum state and the 
oscillator in the first excited state, and states with one photon and the 
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oscillator in the ground state. The survival and transition amplitudes 
relative to more general states factorize in a sum of product of these 
two simpler types. The survival amplitude S(t) has a particularly neat 
form (see Theorem 3.3): 

S(t) = cie-^l*' + cae-i^e-*** + c 3 e" Ae|t| Ei( A e |t | ) . 

In this formula, the complex numbers 

z k := (-l) h u e - ry e , k = 1,2 

coincide with the complex poles of the analytically continued resolvent 
of L+; —A 2 , is the unique negative eigenvalue of L e , and by Ei we mean 
the exponential integral function. The complex numbers Ci, C2 and C3 
depend on the physical parameters. So, in the time evolution of the 
survival amplitude for the bound states of the oscillator it is possible 
to distinguish three different time behaviours. The first term, depends 
on a resonance on the imaginary axis originated from the projection 
onto the positive spectral subspace; it is a pure exponentially decaying 
term (for positive times) and the characteristic time of the decay, for 
realistic values of mass and charge of the electron, has an order of 10 -23 
sec, an exceedingly small time. The second term is an exponentially 
damped oscillation described in terms of the complex resonance poles. 
In particular, the imaginary part of the resonance 7 e gives as usual a 
measure of the lifetime of the excited unperturbed states, or, equiv- 
alently, the breadth of the spectral lines of the spontaneous decay of 
the excited states (according to the Breit Wigner law, see jHj); while 
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the real part gives the position of the maximum in the emission of the 
spectral line. In terms of given physical parameters of the systems, the 
behaviour of these quantities is the following: 

Here the symbols m, c and e denote the (renormalized) phenomeno- 
logical mass, the velocity of light and the electric charge respectively. 
The noteworthy fact is the Lamb shift in the expression of u e . The 
maximum in the emission does not appear in correspondence of the 
unperturbed frequency of the oscillator, but at a displaced frequency. 
The last term, taking into account the asymptotic behaviour of the 
exponential integral, is of the order | for t >> 1. The appearance of 
a slowly decaying tail implies a departure from the purely exponen- 
tial decay given by the Breit Wigner law, and is well known both in 
theoretical models and experimental studies. See Lemma 2.5 and the 
following remarks for an interpretation of its origin. 
The transition amplitudes A(t) have a form very similar to the one 
for S(t) (see Theorem 3.5), and to them apply the same remarks and 
comments concerning their time behaviour. 

2. An Exactly Soluble Model in Classical and Quantum 

Electrodynamics 

2.1. Classical theory. The classical Pauli-Fierz model for a particle 
with charge e, charge density ep r and bare mass m r interacting with 
the electromagnetic field in dipole approximation and subjected to a 
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restoring harmonic force, is described by the hamiltonian 



1 1 



H = 2nc 2 (E, E) + — (A, AA) + - — p - - (p r , A) 



8ir ' 2m r 



e 



c 



2 I 



where E = A/(4irc 2 ) is the canonical variable conjugated to A and p 
is the canonical momentum conjugated to the particle position q. The 
model is written in the Coulomb gauge, so the field are divergenceless. 
A suitable phase space for this dynamical system is H~l(R 3 ) ©L 2 (R 3 ) © 
M 3 ffiIR 3 where Hl(M?) denotes the space of divergenceless locally square 
integrable vector fields with square integrable first derivatives. A point 
electron should have a Dirac measure as charge distribution, and the 
regular form factor p r is introduced to give meaning to the the equations 
of motion. So r has to be interpreted as a measure of the particle radius. 
In the point limit, as r j 0, the charge distribution p r weakly converges 
to So, and the Hamilton equations corresponding to the Pauli-Fierz 
hamiltonian loose their original meaning. A well defined dynamical 
system is recovered only at the expence of renormalizing the bare mass 
m r . This procedure is analyzed in detail in (T^j and we content ourself 
to say here that the correct prescription is given by the well known 
relation between the bare, electromagnetic (m em ) and renormalized 



87r e 2 i 
m = m r + m em = m r + — — - ((-A) p r , p r ) 

3 c z 



Keeping m fixed to the physical value, this is the only choice for the bare 
mass which allows to obtain a nontrivial limit for the Pauli-Fierz model. 
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From now on, with the symbol m r we mean precisely the function 
m r (e, m) given by the above relation. 

The rigorous deduction of the limit dynamics is carried out in [20] 
for a general external potential; specializing the result to the case of 
the harmonic potential one obtains the equations 



A = Anc 2 E 
1 



E = — 

Q = Qa 
p = —aq 



The p-dependent operator is an affine deformation of a linear self- 
adjoint operator in the class of point interactions (see pQ and references 
therein for the use of point interactions in quantum mechanics): is 
linear if and only if p = and is the vector- valued version of one of 
such linear operators; Qa is a certain linear functional which in some 
sense extract the singular part of the vector potential. 

While obtaining a well defined dynamics in the point limit it is an 
interesting and not obvious result, it remains unclear how to quantize 
such a system. So, according to the point of view we adopt in this 
paper, and in view of the Segal quantization of the system, we would 
like to work with (abstract) second order wave equations of the form 



i(t) = -Li{t) , 
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with a suitable self-adjoint operator L. This is not the case of the previ- 
ous system of equation, due to the presence of a tight relation between 
the dynamical variables contained in the definition of the operator H^. 
Going back to the regularized system, the corresponding second order 
equations are not better from this point of view; they are 

1 „• . . Aire „ , . 
—A = AA + Mqp , 

<r c 

e 

m r q = — (p r , A) - aq, 
c 

and the presence of q and A on the right hand side make them not 
of the desired form. So we prefer to give an construction of the limit 
dynamics independent of the one given in [20] • In fact, making resort to 
the hamiltonian regularized equations above, it is simple to overcome 
this problem. Deriving with respect to time the last equation and using 
the third, one obtains that the couple £ = (A,p) satisfies the abstract 
second order wave equation 

where the operator L r e is given by 

Ll(A,p)= (-c 2 AA-—M(p--(p r ,A))p r ,— (p-*( Pr ,A))) 

The analogous result is obtained by means of the canonical transfor- 
mation which exchanges the particle position and momentum (see also 

m 

q= -P , p = Q. 
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The change of dynamical variable from q to p in the second order (or 
lagrangian) formalism is not particularly relevant to the analysis of the 
problems we are interested in, and the interpretation of the results we 
get. For example, the time behaviour of the classical position q, which 
is important in the calculation of the survival amplitude of the bound 
states of the oscillator, is quite simply related to the time behaviour of 
the variable p: up to a constant, the derivative of the momentum gives 
the position, and this relation is preserved in the limit dynamics. 

Another important point to note, is that the operator L£ is a finite 
rank perturbation of the noninteracting operator L = (— c 2 A, 0). This 
simple structure suggests the possibility that the operator U e has a limit 
for r | 0, and being an unbounded operator such a limit, if existing, 
should be sought in the resolvent sense. The calculation of the resolvent 
of the operator L r e is lenghty but elementary. We omit the proof and 
give the result in the following 

Lemma 2.1. For every z £ C±, the resolvent of L r e is given by 



(n e -z 2 r\A, P ) = (gf*A, 



0)-A r (z)R r e (z)(A,p), 



where 



A r (z) 



-1 



Y e ±iz\x\/c 



±\VCIZ > , 




c 2 4n\x 



R r e (z)(A,p) = (Rl r (z)(A,p),R 2 e r (z)(A,p)) , 
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A-7TP 

Rl r (z)(A,p) = M((^ce((-c 2 A-^ 2 )- 1 A,p r ) 

c v 

+c 2 TlK + k r )p \ ( _ C 2 A _ z 2yl p 

a J 
R 2r {z){A,p) = a- ((-c 2 A - z 2 y l A, Pr ) + m r k[p, 



^ = l + -7re 2 ((-c 2 A-^ 2 )-V,Pr) , 

ir a 2 8 a7re .. 2a 2\— i \ 
k 2 = z - -^T7((-c A-z ) p r ,p r ). 

The next step is the point limit r J, 0. The result is analogous to the 
corresponding result given in JH] for the case of a free particle, and 
the proof is modelled on one of the well known ways of defining point 
interactions (see [Tj). We give only an outline of the proof. 

Lemma 2.2. Let p r — ► 5q weakly as r J. 0. For every fixed z G C±, 
(L r e — z 2 )^ 1 converges as r I in the norm resolvent sense to the 
operator 

(gf*A,0)-A ± (z)Rt(z)(A lP ), 
where, putting uj 2 := — ; 

(Aire 
— M(z 2 ce(gT 1 A}+c 2 p)gf, 



±i g z 3 - m(u 2 - z 2 ) 



Such an operator is the resolvent of a self-adjoint operator L e on L 2 (IR 3 )( 
IR 3 when on the component IR 3 one considers the scalar product (pi,p2) ' z 
k-oPi -P2, with k := ^4. 
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Proof. The proof of the convergence of the regularized resolvent is a 
direct consequence of the following limiting relations: 



\im(-c 2 A-z 2 y 1 Pr = g± J 

rj.0 



2e 2 

lim ki = , lim k\m r = m ± z — - z . 

rjQ rlO 3C 3 

limK _ ±i|S^-mK-^) 



Checking that the limit operator is the resolvent of a self-adjoint oper- 
ator is routine. See 1 19 1 for a similar verification. □ 



By the resolvent just constructed it is straightforward to derive the 
actions of L e itself and of its spectral properties: 

Theorem 2.3. The action and domain of the self-adjoint operator 
L e : D(L e ) C Ll(R 3 ) © M 3 -> L 2 (M 3 ) © ]R 3 

are given by 

L e (A,p) = (-c 2 AA , ul (p- ^o(O))) , 

f 47re 
D(L e ) = I (A,p) 6 LffK 3 ) e R 3 : A = A) + — Mi> g , 

V^ , A^L^),„rf,™= P -^ o(( »}, 

w/iere 

0(s) ' 



47r|a;| 
Moreover, 



CTp(L e ) = {-A 2 } , cr ess (L e ) = a ac (L e ) = [0, +oo) , cr sc (L e 
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where 

. 3mc 3 _ o. 

is the unique real (and positive) solution of the third order equation 

^A 3 -mK 2 + A 2 )=0. 

We emphasize that, due to the fact that the operator L e has been 
constructed as a norm resolvent limit of the operator L r e , the flow gen- 
erated by the limit operator coincides with the limit of the regularized 
flow, in the relevant norms of the phase space and uniformly in time. 
This allows to consider L e as the generator of the limit dynamics. As a 
second remark, note that the algebraic equation A-t(z) -1 = is noth- 
ing but the characteristic equation of the Abraham-Lorentz equation 
(see e.g. j!7j ) in the particular case of an harmonic external force, i.e. 

2e 2 



-r Q + q + (J q = r 



3mc 3 



the equation which classically describes the behaviour of the particle 
position in the nonrelativistic regime (the classical relativistic equation 
was obtained by P.A.M.Dirac in JUj). In particular, the negative eigen- 
value of L e corresponds to the so called runaway solution of the A-L 
equation. In the traditional approaches, this solution is discarded due 
to its unphysical character. From the present point of wiew, we give up 
the (important) interpretative problem related to the presence of these 
instabilities, and take the attitude according to which the suppression 
of runaway behaviour corresponds to reduction of the dynamics on the 
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stable subspace, or equivalently, restriction to the absolutely continu- 
ous component of the spectrum. This should correspond, in ordinary 
scattering theory for Schrodinger operators, to the elimination of bound 
states. The procedure to obtain this reduction can be explicitely per- 
formed as follows. Let us consider, from now on, L e as acting on the 
complex Hilbert space L*(IR 3 ; C 3 ) © C 3 . It is self-adjoint when on the 
component C 3 one considers the scalar product (Ci, C2) := «oC* ' C2- 
Given d, ( 2 , C35 an orthonormal base in C 3 , and defining 



be the normalized eigenvectors corresponding to the eigenvalue — Ag. 
Then the projection P e onto the absolutely continuous subspace of the 
operator L e is given by 



We define the positive self-adjoint operator by := P e L e , and 
from now on we consider this reduced operator as the generator of the 
physical limit dynamics. 




let 




i=l,2,3 



In particular 
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2.2. Quantum theory. The operator generates the classical evo- 
lution, whereas, according to the results summarized in the appendix, 
the corresponding quantum evolution is given in terms of its square 
root (L+) 1/2 . More explicitely, denoting by ^(L^M 3 ; C 3 ) © C 3 ) the 
bosonic Fock space over I^(1R 3 ; C 3 ) © C 3 , i.e. 



where S n denotes the symmetrization operator on the n-th sector, the 
quantum hamiltonian on JF(L 2 (IR 3 ; C 3 ) ©C 3 ) corresponding to the sec- 
ond quantization of the classical wave equation 



F(Ll(R 3 ; C 3 ) © C 3 ) := S n (Ll(R 3 ; C 3 ) © C 3 ) 



n>0 



{dl + Lt) (Ap) = 



is given by 



H e := hdT{{L 




Note that the noninteracting hamiltonian 



h := hdT(Ll /2 ) = yr^ewo) 



is unitarily equivalent to 



hdr(V^A) © l + 1 © hdv(uj ) , 



defined on the Hilbert space 



HLliM 3 ; C 3 )) © ^(C 3 ) := S n Ll(R 3 ; C 3 f " © S n (C 3 ) 



n>0 n>0 



The unitary operator 



U : JF(L 2 (R 3 ; C 3 ) © C 3 ) -> JF(L^(IR 3 ; C 3 )) © T(C 3 ) 



(1) 



16 MASSIMO BERTINI, DIEGO NOJA, AND ANDREA POSILICANO 

giving the stated equivalence is denned by 



un-.= n®n, uco)u 



1 



C((p) ® 1 + 1 ® C(C) 



where f2 denotes the vacuum, ip = (<p, () and C is the usual creation op- 
erator. Moreover hdT(uo) is unitarily equivalent to the usual harmonic 
oscillator hamiltonian on L 2 (M 3 ; C) given by the self-adjoint operator 



Concerning the interacting Hamiltonian H e , the present one is its first 
explicit construction in the ultraviolet limit and the problem arises if 
other representations more directly confrontable with usual canonical 
formulation could be given. We emphasize the fact that in our descrip- 
tion, already at the classical level, the point limit produces an intimate 
interlacing between field singularities and particle variables, through 
the definition of the domain of the operator L e itself, and this fact 
introduces essential difficulties in tracing the relation with the canoni- 
cal formalism based on the usual regularized Pauli-Fierz Hamiltonian, 
where this constraint disappears. Nevertheless, in the quoted Arai pa- 
per ([2 ), a reconstruction theorem based on the limit of the Wightman 
functions of the regularized model is outlined. It could be interesting 
to indagate the relations between two approaches. 
Since the resolvent of L e is quite explicit, making use of the Birman- 
Kato invariance principle to deal with the group generated by the 
square root, and of the Birman-Kuroda completeness theorem (see [7]) 



A + 

2m 
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which is applicable because (L e + z)^ 1 — (L + z)~ x is a finite rank 
operator, one immediately obtains the following 

Theorem 2.4. Let P and P e be the orthogonal projections onto the 
absolutely continuous subspaces of L and L e respectively. Then the 
Moller wave operators exist, they are complete and, 

ft ± (H e ,H ) :=s- lim e - ttH «/ R e -* H °/ R r(Po) = r(fi±(L e , L )) 



fi±(H ,H e ) :=s- lim e - itHo/h e- itHs/h T(P e ) = T(tt±(L , L e )) 

A final result of this paragraph is a formula for the evaluation of tran- 
sition amplitudes of the Schrodinger-like propagator e~ l *( Lf ^ 1/2 , that is 
the scalar product of the type (^i, e~^ L ^ 1/2 ip 2 ) , in terms of boundary 
values of the resolvent of the classical operator L e . In the formula the 
special form of the operator L e plays no rule, and it holds true for a 
positive generator A whatsoever. By linearity it will be sufficient to 
suppose that ipi an d ^2 are re &l valued. Since (see the appendix for 
the definition of W L +) 



{{i W Lt _ o) = (z(Lj - z 2 )~V, - iz\Lt - z 2 )' 1 ^) , 
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by Lemma A. 2 and Lemma A. 3 one obtains 



= lim lim / d\e- itx ((^ 1 ,((X + ie)(L e -(X + ie) 2 )- 1 

afoo e|0 2,711 J_ a 

- (A - ie)(L e - (A - ie) 2 ) -1 )^)) 
+ lim lim -L / dA e _iU (Mi, ((A + ie) 2 (L e - (A + ie) 2 )" 1 

atoo ej.0 2ni J_ a 

-(X-ien^-iX-ief)- 1 )^:)-^)). 
Moreover, using first resolvent identity and 

(Ltr i/2 = - 1 dsiLt+sr 1 

n Jr 

it turns that 

(Vi.tie-^r'wr 172 ^) 



Thus for any couple ^i, ^ 2 for which the limits 



(L e - A|)- V 2 + ) := lim (V>i, (L e - (A ± ^e) 2 )" 1 ^ 
40 



exist, one obtains the following 
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Lemma 2.5. 

= lim -L P dX e- itx A ((^, ((L e - A^)- 1 - (L e - A 2 )- V 2 + » 

afoo 27TZ J_ a 

1 f ds f e~ ia A 2 



tc J R 2m J R s 2 + X 2 

Note that the second contribution in the previous formula comes 
from the nonlocal relation between the real phase space classical vari- 
ables and the complexified ones, described in Lemma A. 3. Performing 
however the s integral one obtains, under the same condition of the 
previous result, the following alternative representation 

= lim -L f dX e- itx X ((L e - A 2 ,)" 1 - (L e - Xi)- 1 )^)) 



This last representation has a more direct meaning, in that it is an 
integral extended over the spectrum of the operator (L+) 1 / 2 . Moreover, 
it presents the evolution generated by (L+) 1//2 as a Fourier transform 
of a function (for every fixed couple of states ipi and ip 2 ) supported 
on a half line. As a consequence of the Paley- Wiener theorem, the 
evolution of the amplitude transition cannot have a leading large time 
contribution of exponential type, but a slower one should appear (see 
[El [S] for an early application of this remark to the time decay of 
amplitude transition). The exact time behaviour cannot be precised 
without the knowledge of further details about the generator. To this 
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end it is of more practical use the formula given in Lemma 2.5, as we 
see in the following section. 

3. Radiation theory 

3.1. Generalities. In this section we want to give some details of the 
quantum dynamics of the model we are studying. In particular, we 
want to estimate, under the dynamics generated by H e , the survival 
amplitudes of the bound states of H and the probability amplitudes 
of the transitions between two of such states with emission of photons. 
We begin with some preliminaries on the general structure of the am- 
plitude transitions in our model. 

In the space JF(LJ;(IR 3 ; C 3 )) © J^(C 3 ) a bound state (n-level) for the 
hamiltonian operator Ho is represented by a vector of the kind 

n © s n (Ci ® • • • ® Cn) • 

According to (1) this state is represented in JF(L^(1R 3 ; C 3 ) ©C 3 ) by the 
vector 

S„((0,Ci)<8>---<8>(0,Cn)) 
By (1) again, a more general state with m photons, of the kind 

S m (<Pi © • • • © ip m ) © S n (d © • • • © Cn) , 
is represented in JF(LJ;(R 3 ; C 3 ) © C 3 ) by the vector 



s m+n ((<pi, o) © • • • © (<p m , o) © (o, CO © • • • © (o, • 
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Note that, being H e = hdT((L+) l l 2 ), any sector in JF(L^(R 3 ; C 3 ) ©C 3 ) 
is preserved under e~ lMe . Thus only the transitions from a n-level to a 
m-level with the emission of n — m photons are allowed. 

On the other hand, as we shall see, asymptotically the particle part of 
the wavefunction of the state relaxes to the ground state, and this yelds 
to the asymptotic conservation of the photon number in the scattering 
process. This property was devised by A.Arai in for the regularized 
Pauli-Fierz model with quadratic potential, as a consequence of the 
factorization properties of the scattering matrix, and we find it again 
in the point model. 

Since 

e-^S^i ® • • • ® f/> n ) = r(e-^ 1/2 ) ® • • • ® ipn) 
n! ' 

the survival amplitude of the bound state fi® S n (£i <S> • ■ ■ <2) Cn) is given 
by 

cr J=l 

whereas the the probability amplitude of the transition 
is given by 

II (K,0),e-^) 1/2 (0,0)) II ((0,^),e-^ + ) 1/2 (0,a)). 
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3.2. Survival amplitudes. After these general remarks, we evaluate 
the survival amplitudes for the bound states of the unperturbed dy- 
namics along the perturbed evolution. To simplify the exposition we 
break the analysis in a number of lemmata. 

Lemma 3.1. For any Ci, C2 £ C 3 one has 

((o,Ci),e-^ +)1/2 (o,c 2 ))=5(t)Ci-C 2 , 

where 

s(t) = - 2 -^i(t) 

K 

and, with ±t > 0, 

q(-i\t) z± + ^Ae q'{z±) 



e 






s\ 


p(i\ 


s 


1) 


2(1 


s\ 


+ A e ) q(i\ 


s\ 


)' 



J^t) = - ds 
n Jr 

T(t) = - Id ( eTN<|s| e^X 2 e p(-t\ e ) \ 

2[) vr A S \ 2(\s\ - A e ) q(-i\s\) S 2 - A 2 q(-i\ e ) J ' 



where 



. . e 2 A e . m ( 2 A 2 
. . 2e 2 2 . mcjn , 



z± := ±u e - i% 



e 1 



are £/ie iwo roots ofq(z). 
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Proof. Let z G C±. Taking in account the algebraic equation satisfied 
by A e and the identity 

1 1 



4ttc 3 A e =F iz 
we have, for the projected resolvent, 



((0,(i),(Le-zr\0,( 2 ) + ) 

^2 A S/nT „ s , ( , .2e 2 



= - k — A±(z) \^nuj 2 8^(0$, G Xe ) + \m ± % — zj X e j (* ■ ( 2 

-2^g (±* » * ~ g ("o 2 + f )) + Ae) 
tf(±i*-A.) (g^ + ^ (± ^_ Ae) j (± ^ + Ae) 1 ' 2 " 

Thus using the representation formula given in lemma 2.5, we get 

((o,c 1 ), e -^ + ) 1/2 (o,c 2 )) = -^ (h(t) + - [ dsi 2 (t, S )) cr - c 2 , 

where 

1 /" a 



Jx(t) := — lim / d\ A (/(A) - /(-A)) , 



a 



1 /"^ „-,a^ 2 /(A)-/(-A) 



J 2 (t, s ) := J_ / rfA e- iU A 
2^ Jr 



/(A) := 



s 2 + A 2 ' 

e 2 A e \ I „• m ( , ,2 , *e \ 

i ^j a + 1 a: ^o + tJ _ 1 p(A) 



A + a e |4A 2 + ^(A + 2A e ) A + iAe g(A) ' 



3c3 /v ^ " A| 

Then, by residue theorem and Jordan's lemma, when ±t > 0, one 
obtains 

hit) = T e^iX e P -t^ ± T , 
9(-»A e ) fc 4^ 2 ^ + zA e g'(2fe) ' 
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where Zk :— {—\) k uj e — j e , and 

_ e T\s\t N f 1 P(=P|s|) + 1 P(^kl) 

! \s\ =F A e g(=RN) |s| ± A e g(±«|s|) 



Since 



one has 



7T J R S + £ fc 



+ - I dsl 2 (t,s) 

71" Jr 



e q(-i\ e ) z k + iA e 

+ - I ds 







9(- 





+ - [ ds 





\s\t: 




p(i| 











where fc = 2 if t > and fc = 1 if t < 0. □ 

We distinguish three contributions in the previous formula for I(t). 
The one on the first line is due to the resonances of the system. In 
particular, we have a purely exponential term with lifetime A e , com- 
ing from the projection on the stable subspace of the system, which 
dies out very quickly; and terms exponentially damped, due to the 
complex resonances. These last resonances are typical of the Breit- 
Wigner distribution, and they are the only terms which survive to the 
usual single pole approximation for the resolvent. They give both the 
breadth of emission lines in the spectrum (or equivalently the lifetime 
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of the process) and the Lamb shift, as recalled in the introduction. The 
other contributions to the survival amplitude need more direct analy- 
sis. To simplify the exposition we confine to the case of t > 0, the case 
t < being completely analogous. The term J\{t) is nothing else that 
the Laplace transform (which we denote by the symbol C) a rational 
function. It is a tedious but standard calculation to verify that such 
a Laplace transform is in fact the sum of a pure exponential and a 
damped oscillation with the same characteristic exponents as the ones 
coming from the first group of terms. So, the term J\{t) corresponds 
to a further resonance contribution. 
Concerning J2(t) the following result holds true. 

Lemma 3.2. There exist complex constants Ci,c 2 ,c 3 such that, for 
t > 

m = U ( c ^ ) + «p.v. r sU, 

vr \{s- 7e) 2 + ut) J vr J s - A e 

Proof. Introducing a parameter e to isolate the singularity in A e , one 
can write 

2 [■»*> ( e-^\s\ p{-i\s\) e-^\l p(-iA e ) \ 
2U Wo { 2 (\s\-\ e )q(-t\s\) s* - % q(-i\J J 



2 



e 


M*| 


s\ 


p(-i\ 


s\) 


2(| 


s\ 




A e ) q(-i\ 





= lim [ — [ / ds——. -— — — -—- + / ds 



e 




s\ 


p( 


2(1 


s\ 




Ae) q( 



7T VJo s 2 - A 2 g(-iAe) Ae+6 s 2 - A 2 g(-iAe) 
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and the integrals in the second braces is vanishing for every e. So, one 
remains with 



Finally, the rational function in the above integral can be decomposed 
in a singular (not integrable) part and an integrable one, obtaining 



where Ci, c 2 and c 3 are constants depending on the parameters. □ 

The first integral has the time behaviour of the other already studied 
contributions, and the last one is related to the exponential integral 
function: 



From the well known asymptotic behaviour for large arguments of the 
exponential integral function, one deduce that the leading contribution 



1. In summary, collecting the previous calculations, we end with the 
following result 

Theorem 3.3. There exist complex constants c\,C2,cz depending on 
the physical parameters e, m, c, u such that the survival amplitude is 
given by 






for large times to the survival amplitude is of the order t when t ^> 



S(t) = C\e 




-7e|*l e -^et + C3e -Ae|t| Ei ( Ae |^_ 
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3.3. Transition amplitudes. Now we turn to the evaluation of the 
transition amplitudes. As we saw in the previous theorem, the un- 
perturbed excited states of the oscillator decay. The transition to a 
lower energy state take place with the emission of electromagnetic ra- 
diation, photons in the Fock space. To give a quantitative estimate of 
the probability of emission of photons one has to decide which func- 
tions describe the one particle states of the electromagnetic field one is 
interested in. We choose to calculate the transition amplitudes in the 
case of a regularized plane wave of (approximatively) given momentum. 

Classically, a natural choice for the functions f representing the pho- 
tons of energy fiv should be the divergence-free plane waves of the kind 



Of course such functions are not square integrable and so we consider 
the divergence-free regularization defined by 



(pt(x) := [k — ) A £ e -^Wc-6M/ C) zy > 0, e > 0, |A;| = 1. 



We choose this one as the 1-particle wavefunctions describing photons. 
Eventually, we are interested in removing the cutoff e from the ampli- 
tude transitions. 

We anticipate the definitions of some quantities appearing in the 
statement and proof of the following result. To evaluate the resolvent 
between states relevant to the transitions, one need the scalar product 
(ip e , Qf). Assuming that k is along the z-axis and ( is along the y-axis, 



(p(x) := k A ( e 



—iuk-x/c 



k\ = l. 
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by elementary calculations one obtains 



We have written the function x £ (±z) just defined, as the sum ^(±2) + 
Xi(±z) with two summands, a rational one and a logarithmic one. The 
first, Xri^z) has poles at the points ±a ± ie, and at the same points 
the logaritmic part xK^ 12 ) has branching points. Due to this fact, 
some cautions are needed in the use of the residue theorem to evaluate 
transition amplitudes. Note, moreover, that these branching points of 
X are to be thought as artificial byproducts of the regularization of the 
plane wave, and that the poles of the function Xr correspond to the 
frequencies of the plane wave. These poles depend on the particular 
wavefunction representing the photon, at variance with the poles of the 
function l/q(z), which are determined by the physical parameters of 
the system. We indicate with C the logarithmic cut in the complex 
plane, and we distinguish the components of the cut with real part of 
a fixed sign (±), with C ± . With these premises, we state the following 

Lemma 3.4. For any d, ( 2 £ C 3 one has 



± 




(kAC)(x e r (±z)+xK±z)). 



(o,c 2 )}=-4 e (*)MCi*AC2), 



where 



4-7re K2 



- 2 m 



c 
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and, with ±t > 0, 

F(t) = =fe™ A 2 XHK) ± e^ G ' tUJet2z± 2z± + ' K v(z + ) 
[) T e qHK) z ± + i\ e q'(z ± ) X{Z±) 

eT€t e* ivt 2{±v - ie) 2v + i{X e - 2e) - ie 
±u + i(X e — e) q(±v — ie) 2u 



i 



e TXet x s x (_a e ) eTN'|s|(-2|s| + A e ) x(— 



7r V s2 - A e?(-^e) 2(|s| - A e ) qr(— 

+ - / 



eTl s l*| s |(2|s! 


+ A e ) 




2(| S | 


+ A e ) g(i|s|) 



v z J c z + i\ e q(z) 

+ f dze-** 2z + iK — - f dze-^ 2z + iK —) 

z/ 3 \J C+ z + i\ e q(z) Jc- z + i\ e q(z) J 

Proof. One has, proceeding as in the proof of Lemma 3.1 



(( lf l,0),(L e -z 2 )- 1 (0,C2) + ) 

Alices (±2iz — \ e )(±iz + A e ) 



x \±z) (*Acr)-c2. 



K*&Z + A e ) (g ^ _ ^ (:pi;8 + Ae) ) (± fc + Ae) 

Therefore 

((^,0),e-^ + ) 1/2 (0,C 2 )) = -— ^ + - / da/afts)) (*A£K 2 , 

where (we omit the dependence on e) 

h(t) := -L lim f dA e"^ A A ((?(A) - «?(-A)) , 

2m a|oo ,/_ 

J 2 (M):= — [ d\ \> 9{X) ~ 9{ ~ X) 
Z^i J R 



s 2 + A 2 
2A + iA e x(A) 



A + iA e g(A) 

Now, let us choose a path in the complex lower halfplane wich has the 
real axis as the upper side, avoids the cuts of the function xi ( sa y> the 
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straight half lines parametrized as z = ±u + u — ie , ±u > 0), and close 
itself at oo along a great circle, as in the previous lemma. There will be 
contributions due to the residues of the function l/q(z), the residues 
the function Xr( z ), an d to the discontinuity of the logaritmic part of 
Xi( z ) along the cut. In the end, one obtains, when ±t > 0, 



9(-*A e ) *-f 2 z k + i\ e q'(z k ) 



where z± are the poles of q, w k '■— {—\) k u — ie are the two poles of 
Xr( z ) and r k the residues of Xr, and finally C is the path along the cut 
of the logarithmic term. 
Moreover we have 



hit, s) = i 



S 2 _ A 2 q{ _ tXe) 



+ e 



=F7e< 



E 



z k 



-f 2 s 2 + zl z k + i\ e q'(z k ) 

Tet w k e Tiut w k 2w k + i\ e 

+ e 1^ ~2~, — 2 — ? — r~ r ' k 

* i 2 + w k w k + «A e q{w k ) 



X( z k) 
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s| 


+ A e 
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|^A e g(=Fi|s|) 



,2^ + iA e ^ 2 



z + iA e 5(2;) s 2 + z 2 
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Therefore one obtains 
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hit) + - I dsl 2 (t,s) 



71 



q{-i\e) z k + i\ e q'{z k ) 

Tet e^2w k 2w k + t\ e 
w k + i\ e q{w k ) 

+ - I r/.s I — — f ^ 4 



e =F|s|* 




-2! 


s 


+ A e ) 


2(| S | 




Ae) 



7T Jr \ s 2 - \ 2 e q(-i\ e ) 2(|s|-A e ) q(-i\s\) 



+ - ds 
71 Jr 



e 

+ 



e=FN*| s |(2|s| + A e ) 




1*1) 


2(|s| + A e ) 


q(i\ 


*\) 


itz 2z + i\ e 1 







v 3 J c z + i\ e q(z) 

e ( f , -itz 2z + iX e Z f Uz 2z + i\ e Z 



+ TT— dze- ltz / dze 



v 3 \Jc+ z + i\ e q(z) J c - z + i\ e q(z) 

where k = 2 if t > and k — 1 if t < and are the components of 
C with ±Sign(Re(z)) > respectively. 

□ 

Now we can give the time behaviour of the various terms, as in the 
previous theorem. A first group of terms is composed by the resonant 
contributions. One distinguish natural resonances, depending on the 
structural parameters only(the physical constants e, m, c), and the 
resonance due to the incident photon. For e J, this last contribution 
reduces to a strictly oscillating term, as expected. A second group of 
terms is given by the s integrals. We give their behaviour for vanishing 
e only. The calculation is similar to the one already given for the sur- 
vival amplitude, and one has in the end Laplace transforms of rational 
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functions with poles at the resonances, producing other exponentials 
of the type already seen, and an exponential integral function. So the 
leading behaviour of the type t~ l survives to the e j limit. 
The last group of terms are the contour complex integrals. These could 
be analized asymptotically as Fourier integrals of rational functions, 
but they vanish as e j 0. Summarizing we can state the following 

Theorem 3.5. There exist complex constants C\, C^, C3 depending on 
the physical parameters m, e, c, u , v , and a function R(t) such that 

lim A e (t) = Cie^ 1 ' 1 + Cae-^'e-*** + C 3 e~ iut + R(t) 

with 

R(t) = 0(l/t), \t\ >1. 

Appendix A. Quantization of abstract wave equations 

In this appendix we give a brief and self contained introduction to 
the quantization of second order abstract wave equations, along the 
lines traced by I.E.Segal in the fifties and sixties and with an emphasis 
on the aspects of direct concern with our work. We refer for details and 
different approaches to [27] and [T5j as regards abstract wave equations 
and to [22], as regards quantization. 

Let B : D(B) C 7i — > TC be a injective self-adjoint operator on the 
real Hilbert space (Tl, (■, ■)). We denote by Tt 1 , the Hilbert space given 
by the domain of B with the scalar product (■, - )i leading to the graph 
norm, i.e. 

(01, <p2)l ■= (Bfa, B<f) 2 ) + <j) 2 ) . 
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We then define the Hilbert space 7Y 1 by completing the pre-Hilbert 
space D(B) endowed with the the scalar product 

[0i,0 2 ]i: = (BfaBfa). 

We extend the self-adjoint operator B to H 1 by considering B : H 1 — > 
Tt, the closed bounded extension of the densely defined linear operator 

B : H 1 C H 1 -> W. 

Since -B is self-adjoint one has Ran(_B)- L =Ker(_B), so that, being _B 
injective, Ran(_B) is dense in H. Therefore we can define B' 1 : H — > 
as the closed bounded extension of the densely defined linear operator 

B- 1 : Ran(B) CH^H 1 . 

One can then verify that 5 is boundedly invertible with inverse given 
by B-\ 

Given B we can now introduce the space 7Y 2 defined by 
H 2 := {0 e U x : B(f) G H 1 } . 
On the Hilbert space H 1 ®H with scalar product 

(((01, 0l), (02, 02))) := Wl, 50 2 ) + (0 1; 2 ) . 

we define the linear operator 

Wb-.^QH 1 CH 2 QH^H 1 ®H, W b {<P, 0) := (0, -BBcj)) . 
Theorem A.l. The linear operator Wb is skew-adjoint. 
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Being Wb skew-adjoint, by the Stone theorem it generates the strongly 
continuous one-parameter group of isometric operators e tWB . By defin- 
ing (<f>(t), 0(t)) := e tWB ((f), 0), one has that, in the case G H 2 , G H 1 , 
4>{t) is the unique strong solution of the Cauchy problem 

'(pit) = -BB(f)(t) 
0(0) =0, 0(0) = 0. 

Given an arbitrary real Hilbert space (/C, (•, -)k:) we will denote by 
(/C c , (•, -) c ) its standard complexification, i.e. /C c = /C © /C, the multi- 
plication by the complex unity being defined by iip := Jip, J(y?i, (fi 2 ) : = 
(-<f 2 ,<fi), and 

Given any linear operator A : D(A) C /C — > /C on /C, we define 
A c : -D(A C ) C /C c -> /C c by D(A C ) := D{A) x .0(A), A c (ipi, <p 2 ) := 
(A<fi, A<f 2 )- Conversely, given any linear operator L : D(L) C /C c — > 
/C c , L((p U (p 2 ) = (Li(v?i, v? 2 ),^ 2 (v ? i> ^2)), we define L r : D{L r ) C /C -> 
/C by D(L r ) := {0 G K : (0,0) G L r := Li(</?,0). 

By the above definitions _B C is self-adjoint and (Wb) c = Wb c , where 



By the Stone formula one has the following 
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Lemma A. 2. 

e tw B 

= limlim / dX e- itX ((iW Bc - (X + te))' 1 - (iW Bc - (A-^e))- 1 ),. 

afoo e|0 2m J_ a 

The following lemma translates the wave flow in a Schrodinger like one. 
Lemma A. 3. The map 

Cb-.H 1 ®H^H c , C b (0,0) := (50,0). 

is unitary once one makes Ti 1 © Ti a complex Hilbert space by intro- 
ducing the complex structure Jb{4>, 0) : = {—B^ 1 ^,!^^) and by defining 
i(0, 0) '■— Jb(4>,4>) and the scalar product 

(((0!, 0x), (0 2 , 2 )» B := (((01, 0l), (02, 02))) -i(((01, 0l), ^(02, 2 )» ■ 

Moreover 

e w s = c * b e -itB c Cb _ 

Once we have transformed the abstract wave equation <p(t) = —BB<f)(t), 
defined on the real Hilbert space H 1 (BTC, into the Schrodinger equation 
iip = B c ip, defined on the complex Hilbert space 7i c , we can then (sec- 
ond) quantize it in the standard way. Let us define K to be the bosonic 
Fock space over 7i c (see j22], section II. 4, for the definition). For any 
ip G TC we define the self-adjoint operator on K given by the Segal field 
S(^) := ^(C(-0) + C*(?/>)), were C and C* denote the usual creation and 
destruction operators (see [23], section X.7, for the definition). Given 
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the Segal field S we can then define the Weyl system W(^) := e l5 ip, so 
that 

W(?/>i + if> 2 ) = ^ ( ^'^ >c W(V'i)W(^) • 

The unitary strongly continuous one-parameter group of evolution on 
K defined by U(£) := T(e~ ltBc ) satisfies the relations 

\N{e- itBc tp) = U(t)W(V>)U(t)* 

and we denote by dT(B c ) the self-adjoint operator on K with generates 
U(t) (we refer to [2H|, section II. 4, and jU] for the definitions of T and 
dT). The quantum hamiltonian corresponding to the (second) quan- 
tization of the abstract wave equation <p(t) = —BBcj)(t) is defined by 
H := hdT(B c ). Suppose now that we start with a selfadjoint operator 
A on the real Hilbert space /C. Denoting by A + the positive part of 
A and by /C + the projection of /C onto the positive spectral subspace. 
Then we can apply the previous construction to B := (A + Y^ 2 , con- 
sidered as an injective self-adjoint operator on Tt := fC + . However, 
since 

1 /2 

where ip + denotes the projection of on /C + , we will work with dT{{A + ) c ) 
on the bosonic Fock space over /C. 
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